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ABSTRACT 


•i, 


This  research  memorandum  describes  a  queueing  model 
used  to  analyze  a  sparing  decision  for  a  part  with  general 
probabilistic  demand.  It  also  describes  an  extension  of  the 
model  to  include  discriminating  treatment  of  the  repair  and 
resupply  pipelines.  The  final  section  applies  the  model  to  an 
illustrative  example. 


\ 


INTRODUCTION 


Occasionally,  when  a  sea-based  aircraft  is  being  repaired  or  inspected,  a  certain 
part  needs  to  be  replaced  A  prevailing  assumption  1  and  '2  is  that  the  number  of  such 
occurrences  during  a  specific  time  period  can  be  approximated  with  a  Poisson  probability 
distribution.  An  underlying  assumption  leading  to  this  probability  is  that  the  expected 
number  of  failed  parts  depends  only  on  the  number  of  flight  hours  in  that  time  period 
Recent  research  at  the  Center  for  Naval  Analyses  (CNA)  suggests  that  this  may  not  be 
the  case  3.-1.  and  -5 


Another  characteristic  of  the  Poisson  probability  distribution  is  that  there  is  a  posi¬ 
tive  probability  that  a  very  large  (theoretically:  infinite)  number  of  parts  will  fail  during 
a  large  enough  time  unit.  In  practice,  however,  the  probability  that  the  number  of  failed 
parts  will  be  greater  than  the  total  number  of  parts  (used  or  spared)  equals  zero. 


This  paper  provides  a  tool  for  analyzing  a  sparing  policy  for  a  part  for  which  the 
demand  is  not  necessarily  Poisson  and  for  which  there  is  a  limit  on  the  number  of  parts. 
This  tool  can  be  used  as  the  single-item,  single-echelon  statistical  analysis  component 
of  a  more  complex  optimization  model.  The  next  section  formulates  a  queueing  model 
with  exponential  service,  ample  servers,  and  finite  capacity  The  model  s  general  appli¬ 
cation.  extension  to  hyperexponential  service,  and  application  to  a  numerical  example 
are  considered  in  the  subsequent  section,  and  the  last  section  contains  some  concluding 
remarks 
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MODEL 


When  an  aircraft  is  being  repaired  or  inspected  and  a  ertam  part  is  found  to  have 
failed,  it  is  replaced  immediate!) .  if  a  spare  pari  is  available  Tl.t  availability  of  spares 
is  determined  by  (1)  the  number  of  spares  initially  ■storked  or,  -dot:  ,  Ji  the  nuniber  of 
"new"  spares,  that  is.  the  number  of  parts  repaired  or:  ship  or  ordered  ana  received 
from  shore:  and  (3)  the  demand  for  spares  already  satisfied  that  is.  the  number  of  parts 
already  used  to  replace  failed  parts  Simultaneous  with  the  rise  of  the  demand  for  spare', 
the  removed  part  enters  a  queueing  system  (figure  1) 

In  queueing  theory  terminology  6  and  7  .  the  removed  parts  are  customers  and  they 
arrive  to  receive  service  provided  by  serrers.  A  five-part  descriptor.  A  B  C  F  A'.  is 
used  to  summarize  symbolically  the  information  regarding  a  queueing  system:  A  is  the 
interarrival-time  probability  distribution.  B  is  the  service-time  probability  distribution. 
C  is  the  number  of  servers.  F  is  the  system’s  capacity,  and  A  is  the  size  of  the  customer 
population.  .4  and  B  take  on  values  from,  among  others,  the  following  set  of  symbols 
w  hose  interpretation  is  given  in  terms  of  distributions  w  ithin  parentheses:  XI  (exponen¬ 
tial).  En  (n-stage  Erlangian).  H_\  (N-stage  hyperexponential).  D  (deterministic),  and  G 
(general).  C.  F .  and  K  are  either  integers  or  infinity,  with  infinity  (oc)  symbolizing  the 
fact  that  whatever  the  actual  number  is.  it  is  large  enough  to  yield  a  certain  attribute 
of  the  system 

In  this  notation,  the  currently  used  model — the  single-item,  single-echelon  compo¬ 
nent  of  the  Multi-Item  Multi-Echelon  (MIME)  inventory  model — is  an  M  XI  x  x  x 
queueing  system  That  is.  the  interarrival  and  service  times  are  exponential:  the  number 
of  servers  is  large  enough  so  that  no  customer  has  to  wait  for  service  to  begin,  the  system 
is  uncapacitated:  and  the  size  of  the  calling  population  is  iarge  enough  for  the  demand 
rate  to  remain  constant,  regardless  of  the  number  of  customers  already  in  the  system 
The  number  of  customers,  i  e..  failed  parts,  present  in  this  system  has  a  non-zero  proba¬ 
bility  of  being  larger  than  the  total  number  of  parts  To  remedy  this  situation,  one  can 
use  the  XI  XI  x  m  oc  queueing  system,  with  m  denoting  the  total  number  of  parts 

It  is  of  interest  to  note  that  the  finit e-number-of-parts  issue  was  addressed  by  a 
finite-system-capacity  model  rather  than  by  a  finite-calling-population-size  model  The 
reason  for  this  action  is  that  the  model  needed  to  reflect  the  fact  that  the  demand  rate 
remains  constant .  regardless  of  t  he  number  of  customers  already  in  t  he  sy  stem  Although 
the  main  reason  for  a  constant  demand  rate  is  that  operable  aircraft  must  also  complete 
missions  of  inoperable  aircraft,  it  is.  perhaps,  enhanced  by  cannibalization  —  an  issue  not 
explicitly  addressed  in  this  paper 

Although  the  XI  XI  oc  rn  oc  finite-capacity  system  is  an  improvement  over  the 
XI  XI  x  oc  oc  uncapacit ated  system,  it  still  lacks  the  generality  to  accommodate  any 
interarrival-t ime  distribution  Thus,  the  G  XI  x  m  x  queueing  system  is  considered 
next  Because  no  more  than  m  customers  can  be  present  in  the  system,  there  is  no 
loss  of  generality  in  assuming  that  only  n  servers  provide  service  as  long  as  n  assumes 
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a  value  from  the  set  { m.m  —  1 . oc}.  In  particular,  when  n  =  m.  the  result  is  the 

G  M  m  m  oc  queueing  system. 

This  system,  along  w  ith  its  special  cases,  is  well  solved  for  the  exponential  (G  =  M ) 
case  6  and  7  .  The  only  explicit  attack  on  the  G  \t  m  m  oc  system  is  by  Takacs 
(see  the  chapter  on  Telephone  Traffic  Processes  of  his  book  S  ).  Although  his  result  is 
significant  in  the  theoretical  sense,  its  use  is  limited  b\  the  fact  that  it  is  given  as  the 
Laplace  transforms  of  the  probabilities  of  interest  That  is.  to  obtain  the  probabilities, 
the  Lapiace  transforms  must  be  inverted  This  cannot  be  done  in  a  mechanized  manner 
The  work  described  in  this  paper  follows  the  solution  approach  to  the  infinite-capacity, 
n-server  case  6  and  7  .  Results  are  expressed,  as  a  matter  of  convenience,  in  terms  of 
Laplace  transforms  of  the  interarrival  time  density,  but  they  do  not  require  inversion. 
The  appendix  shows  that  the  solution  here  reduces  in  the  exponential  case  to  the  well- 
known  results.  The  remainder  of  this  section  is  devoted  to  a  solution  procedure  for  the 
G  M  m  m  oc  queueing  system. 

Because  the  exponential  service  time  is  memoryless,  at  each  arrival  instance  the 
system  starts  from  the  beginning,  in  the  sense  of  statistical  behavior.  Thus,  the  arrival 
instances  constitute  the  regeneration  points  of  an  imbedded  Markov  chain.  The  states 
of  this  chain  are  the  number  of  customers  found  in  the  system,  and  the  probabilitv  of 
transition  from  state  i  to  state  j  is  the  conditional  probability  of  the  n  -  lft  arrival  to  find 
j  customers  in  the  system,  given  that  the  n:h  arrival  found  i  customers  in  the  system. 
If  qr.  is  the  number  of  customers  found  in  the  system  by  the  n:n  arrival,  then  {<?.,}  is  a 
discrete-state  Markov  chain  and  its  transition  probability  is  defined  by  equation  1  as 

Pi]  =  Pr{qn-\  =  j  <?„  =  |}  (1) 

This  transition  probability  is  nothing  but: 

p,:  -  Pr  { ( i '  -  j)  customers  are  served  during  an  interarrival  time}  (2) 

Pr  { ( i'  -  j)  customers  are  served  during  an  interarrival 
time  of  length  t}  dA(t)  . 

which  equals 


j;,x  (’')  (1  -  f  ul  e  -*’-dA(t):i  =  0.  1 


.  m:  j  —  0.1.  .  .  /' 

otherwise  . 


where  i'  =  min{i  -  l.rr?l  and  the  integration  is  Riemann-5t leltjes  It  should  be  noted 
that  p =  p„.  _  ]  ;  for  all  j 
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Figure  2  present?  a  diagram  with  the  transition  probabilities  into  and  out  of  state  t. 


Replacing  ( 1  -  e  )'  J  with  its  hinomial  expansion  ^2'k~,J  |  '  j  (  -  I  )*(e '  u:)“. 

and  not  ins  that  [t~  dA(t)  is  the  expected  value  of  with  respect  to 

random  variable  t  (that  is.  t  he  Laplace  transform  of  d . 4i  I !  evaluated  at  p(  A  -  j ) .  .4 '  p( A  ~ 
f  j  ]  j.  equation  3  can  be  rewritten  m  the  follow  ms.  computationally  tractable,  form 


=  I  (J')  (k  -  j)n  :  i  =  0.1 . m:j  =  "  1 . « '  (4) 

|  0  :  otherwise  . 

w  here,  as  before,  r1  ~  minji  —  1.  rn} 

Let  rk  -  Pr{qn  -  A-}  denote  the  (steady-state,  unconditional)  prcbabilitv  o>f 

a  random  arrival  to  find  A  customers  in  the  system  It  should  be  noted  that  finite 
system  rapacity  ensures  that  steady-state  equilibrium  is  attainable.  Having  found  the 
(one  step,  conditional)  transition  probabilities  given  by  equation  4.  rk  can  be  found  by 
solving  the  following  (m  -  1)  by  (m  -  1)  linear-equations  system: 


rn 

r *  -  V  rnpnk  :  A  1.2 . m 

r.  -  K  1 
m 

Vrn  =  1 

n  -  ■ < 


(6) 


Toward  this  end.  equation  5  can  he  rewritten  as. 


rk 


*ry 


ri  :  K 


Pk  1  K 


A  1.2 . rn 


If.  based  on  observations,  theoretical  knowledge  or  computational  instability  r,„ 
equals  0  (or  is  small  enough),  rn  can  be  replaced  by  rn  -  1:  thus,  there  is  no  loss  of 
generality  in  assuming  r„  •  0.  and  the  following  ratio  rati  be  defined: 


A-  implied  by  equation  6.  r„.  must  equal  1  .  cK:  therefore  equation  8  can  be  rewrit¬ 

ten  as 


n 


That  is.  given  t he  transition  probabilities  p,  calculated  in  the  previous  step,  the  system 
of  linear  equations  can  be  solved  by  finding  ck:  k  —  0.  1..  . .  m.  This  is  done  by  the 
following  recursive  formulae 


=  1 


Ck 


-  T.  c*p 


m  -  n  m  -  k 


(10) 


Ck~\ 


Pm  -  f  4  -  1 )  m-k 


APPLICATION 


An  application  of  the  model  developed  in  the  previous  section  must  address  the 
following  issues: 

•  What  are  the  queueing  system's  characterist ics.  or  what  can  be  said  about  service0 

•  What  is  the  input  to  the  queueing  system,  or  what  can  be  said  about  arrivals? 

•  What  is  the  significance  of  the  queueing  system's  output  to  decision  makers0 


SERVICE 

Because  the  model  presented  here  uses  queueing  theory,  the  information  about  the 
number  of  customers  served  during  a  specific  time  period  is  conveyed  as  information 
about  the  time  it  takes  to  serve  a  customer  Because  a  customer  is  a  removed  part  and 
service  includes  off-ship  resupply  (when  on-ship  repair  is  not  possible),  it  is  natural  to 
consider  service  time  in  calendar  time  rather  than,  for  example,  flight  hours  or  number 
of  flights.  Average  service  time  equals  average  on-ship  repair  time  weighted  by  the 
probability  of  repair  aboard  ship,  plus  the  average  off-ship  resupply  time  weighted  by  the 
probability  of  off-ship  resupply.  For  example,  in  the  case  of  a  pulse  decoder,  5  percent  of 
the  removed  parts  can  be  repaired  aboard  ship,  but  the  remaining  95  percent  must  be  sent 
off  ship  The  average  on-ship  repair  time  is  4  days,  while  the  average  off-ship  resupply 
time  is  44  days.  Thus,  the  average  service  time  equals  42  days  (  05  >  4  -  .95  >  44  =  42). 
This  holds  true  for  any  probability  distribution  of  service  time.  To  facilitate  numerical 
comparison  with  the  currently  used  mode!  and  for  the  sake  of  ease  of  exposition,  service 
time  is  assumed  to  have  an  exponential  distribution  (with  parameter  /u  that  equals  the 
reciprocal  of  the  average  service  time)  Such  a  combination  of  on-ship  repair  and  off-ship 
resupply  implicitly  assumes  that  a  part  repaired  on  ship  is  identical  to  a  shore-supplied 
part . 


ARRIVALS 

The  information  pertaining  to  the  number  of  arrivals  is  conveyed  again  as  informa¬ 
tion  about  interarrival  times  The  choice  of  interarrival  times  to  be  in  calendar  time  is 
not  a  natural  one.  because  it  necessitates  “stretching”  flights  into  calendar  time:  how¬ 
ever.  that  choice  is  dictated  by  the  need  to  be  compatible  with  service  times  and  to 
allow  consideration  of  factors  contributing  to  the  rise  of  demand  for  spare  parts — that 
is.  arrivals — that  do  not  depend  only  on  flight  hours.  Two  approaches  to  determining 
the  probability  distribution  of  interarrival  times  are  possible,  parametric  and  nonpara- 
metric  In  the  parametric  case,  a  probability  distribution  is  determined  by  theoretical 
knowledge  of  the  underlying  processess  or  by  closeness  of  shape  to  the  shape  of  real-life 
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data,  the  parameters  are  estimated  using  reai-iife  data.  In  this  case,  the  Laplace  trans¬ 
form  of  the  appropriate  (theoretical)  probability  density  function,  with  the  appropriate 
parameters  is  used  in  equation  4.  In  the  nonparametric  rase,  an  explicit  consideration  of 
a  probability  distribution  and  its  parameters  is  by-passed,  and  the  data  are  used  directly 
to  calculate  the  Laplace  transform  in  the  following  way  Assume  that  n  observations  of 
interarrival  (calendar)  times  are  available,  from  either  real  life  or  a  simulation.  Denote 
these  observations  by  tj.t*, . . .  ,t„  The  Laplace  transform  in  this  case  is: 


-4-  (*-7V  . 

n 


If  n  is  too  large  a  number,  observations  can  be  grouped  together,  creating  an  empirical 
distribution,  and  in  equation  11  replacing  t,  by  the  middle  of  the  group's  range,  n  by 
the  number  of  groups,  and  £  by  the  proportion  of  observations  falling  in  that  group.  (It 
should  be  noted  that  information  can  be  lost  while  grouping.) 


OUTPUT 

Probabilistic  knowledge  of  the  number  of  parts  found  in  the  queueing  system  allows 
for  probabilistic  knowledge  of  the  number  of  parts  available  for  operational  use.  In 
particular,  the  rks  can  be  used  to  calculate,  for  a  given  outfitting  level  s.  the  expected 
number  of  back  orders  found  at  arrival  points  B(s).  and  the  probability  that  demand 
can  be  fulfilled  immediately  upon  arrival,  i.e..  the  fill  rate  F(s): 


J2  (*  -  s)r* 


F(0)  =  0 


Before  proceeding  to  the  numerical  example,  let  us  revisit  the  assumption  that  service 
time — which  includes  on-ship  repair  and  off-ship  resupply — has  an  exponential  distri¬ 
bution  For  the  currently  used  model,  this  issue  is  purely  theoretical,  its  results  hold 
regardless  of  the  actual  distribution  of  service  time.  For  the  model  proposed  here,  the 
matter  is  crucial  only  because  of  the  memoryiessness  of  the  exponential  distribution 
were  we  able  to  find  regeneration  points,  an  embedded  Markov  chain,  and  to  complete 
the  analysis.  In  a  queueing  system  with  two  types  of  servers,  each  type  having  an  expo¬ 
nential  service  time  (with  different  parameters),  the  combined  service  time  is  exponential 
only  w  hen  all  servers  are  busy.  In  this  case,  there  are  two  Poisson  outputs  of  completed 
services,  whose  combination  is  a  Poisson  output.  In  this  model,  in  general,  not  all  servers 
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are  busy  Assuming  that  each  type  of  service  (on-ship  and  off-ship)  can  be  approximated 
as  having  an  exponential  distribution,  the  combined  service  can  be  viewed  as  having  a 
two-stage  hyperexponential  distribution.  This  approach  is  similar  to  the  one  taken  in 
9  . 

Let  p'  be  the  probability  of  repair  aboard  ship,  p'  is  estimated  by  p.  w  hich  equals  the 
number  of  parts  removed  that  can  be  repaired  aboard  ship,  divided  by  the  sum  of  the 
number  of  removed  parts  that  can  be  repaired  aboard  ship  and  the  number  of  removed 
parts  that  cannot  be  repaired  aboard  ship.  Let  pl  denote  the  reciprocal  of  the  average 
on  -ship  repair  time,  and  p?  denote  the  reciprocal  of  the  average  off-ship  resupply  time. 
The  density  function  of  the  combined  service  time  is  given  by: 


/12(f)  =  PV\e  —  ( 1  ~  p)H'.t  M:'  •  (14) 

Equation  14  can  be  interpreted  as  a  conditional  density  Thus,  the  resulting  G  H?  oc  m  oc 
queueing  system  can  be  analyzed  in  the  following  way: 

1.  Replace  p  of  equation  with  pj . 

2.  In  equations  a  through  10.  replace  rk  with  r\. 

3.  Repeat  steps  1  and  2.  this  time  with  p;  and  rjj 

4  Let  rk  =  pr£  -  (1  -  p)r\.  k  -  0.  1 . m. 

To  illustrate  the  subject,  consider  the  following  numerical  example  There  are  four 
aircraft,  each  using  one  unit  of  the  part  of  interest.  Each  aircraft  flies  once  a  day  for 
four  hours  The  observed  average  failure  rate  is  two  units  per  100  flight  hours.  With  an 
average  flight  length  of  four  hours,  this  translates  into  eight  failures  per  100  failures  are  6 
to  3  to  1.  respectively,  for  equal  times  The  proportions  of  time  spent  in  each  environment 
are  0.5  24.  19.5  24.  and  4  24.  respectively.  The  observed  average  service  time  is  five 
days  for  on-ship  repair,  which  happens  20  percent  of  the  time.  In  the  remaining  SO 
percent  of  the  time  when  the  repair  cannot  be  done  on  ship,  the  average  resupply  time 
is  20  days 

Four  ways  to  model  demand  for  spares  are  considered  here: 


•  An  infinite  number  of  parts,  exponential  interarrival  times,  denoted  I.N’EX  This  is 
the  currently  used  model  and  is  considered  in  order  to  facilitate  comparison. 

•  A  finite  number  of  parts.  (4  ■  1-2  =  6).  exponential  interarrival  times,  denoted 

FINE 

•  Demand  derived  from  four  flights  per  day.  8  100  probability  of  failure  per  flight, 
denoted  LEVY  3  . 


•  A  finite  number  of  part?  and  hyperexponent lai  lriterarrival  times  (a  mixture  of 
three  exponentials,  each  representing  one  environment  of  failure),  denoted  HYPE 

Two  ways  to  model  service  times  are  considered: 

•  Exponential  with  parameter  (0  2  ■  5  —  0  8  -  20)"1  =  IT-1,  denoted  by  EXPO 

•  Hyperexponential  (a  mixture  of  two  exponentials,  one  with  parameter  o~U  the 
other  with  20 —  1 ) .  denoted  HYEX 

The  results  for  all  eight  cases  are  summarized  in  terms  of  expected  back  orders  and 
fill  rates  in  table  1.  It  should  be  noted  that  in  spite  of  mean  rates  being  equal,  different 
distributions  (for  demand  or  service)  yield  different  results.  (For  INEX.  the  service 
distribution  does  not  affect  the  expected  back  orders  )  The  infinite-number-of-parts  case 
yields  the  largest  number  of  back  orders  and  the  smallest  fill  rates.  This  is  not  surprising, 
since  "parts"  keep  failing  after  all  real  parts  have  failed. 

The  closeness  between  FIXE  and  LEVY  can  be  explained  by  recalling  the  Poisson 
approximation  to  the  Binomial.  The  better  results  produced  by  HYPE  are  somewhat 
surprising,  and  they  hint  to  the  value  of  having  a  better  understanding  of  the  failure 
process. 

table  i 

NUMERICAL  RESULTS 
Demand 


Service 

INEX  FINE 

LEVY 

HYPE 

Expected  back  orders 

EXPO 

HYEX 

3.472  2.265 

3  472  2.108 

2  263 

2.104 

1.986 

1.878 

Fill  rate 

EXPO 

HYEX 

0279  .0403 

.1148  .1233 

.0384 

1241 

.0788 

.1549 

11 


CONCLUSION 


In  summary,  the  model  presented  in  this  paper  can  analyze  the  consequence  of  an 
initial  outfitting  level  for  a  part  with  general  probabilistic  demand  a  two-step  hyperex- 
ponential  service,  and  a  limit  on  the  number  of  parts.  The  two-stage  service  is  readily 
extendable  to  any  integer-stage  hy  perexponential  service. 

In  terms  of  continuing  research,  a  critical  area  is  the  programming  of  the  method 
presented  in  this  paper  and  its  incorporation  into  a  multi-item  multi-echelon  inventory 
model  These  areas  will  tie  reported  on  in  a  future  research  memorandum. 
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APPENDIX 


THE  EXPONENTIAL  CASE 


When  the  distribution  of  the  interarrival  times  is  exponential  with  parameter  A.  A{t) 
equals  1  -  < " Al  for  t  0  and  vanishes  otherwise:  the  Laplace  transform.  .4'.  evaluated 
at  ( k  -  j ) is  given  b\ : 


-4-  (A--7> 


A 

A  -  (k-  j)n 


( A — 1 ) 


Dividing  the  numerator  and  denominator  of  the  right  hand  side  (RHS)  of  equation 
A  1  by  p  and  denoting  A  p  by  p  yields: 


A  (A--J)p 


P  _  P  1 

p  -  k  -  j  p  -  J  JL-  -  1 


(A-2) 


Substituting  the  RH?  of  equation  A-2  for  A  [k  -  j)p  in  equation  4  of  the  main  text 
\ields  the  following  transition  probabilities  for  the  exponential  case: 


P<;  = 


J 


-  J 
k 


(-1)* 


_P_ 
P~  J 


1 


(A-3) 


Equation  A-3  (and  its  successors  in  the  sequel )  is  valid  for  t  —  0.1 . m:  j  -  0.1. 

where  i'  -  min{i  -  l.m}.  and  vanishes  otherwise. 

Now.  let  us  concentrate  on  the  rightmost  part  of  equation  A-3 


Replacing  the  rightmost  component  of  equation  A-3  with  the  RH5  of  equation  A-4 
and  moving  -jy  to  the  left  of  the  summation  yields: 


( •: ) 

1  p  V  ( 

i  j  ) 

o  -  J  ^  ' 

r 

1  *■'  ] 

1  e  r  1 

'  j  i 

_ 1 

—> 

1 

J  '(-!)*  (x-) 


k= 0 


di 


(A-5J 


/  1  \  *  ~~  J 

The  bracketed  expression  is  the  binomial  expansion  of  -  1)  .so  equation 

A-5  reduces  to: 


{ ']  )  rn  f' (”L‘ - ‘di  ■ 


(A-6) 


tpon  substitution  of  ( 1  -  i for  t.  and  (p  -  j)(  1  -  y),'~-  ~  Xdy  for  di .  equation  A-6 
yields 


p,j  -  (  ’  )  — —  [  y'  ](p  -  j){  1  -  yV  J~]dy  ( A-7 ) 

\  J  I  P  -  J 

Replacing  y  with  -  ( 1  -  ;)  and  dy  with  d:.  noting  that  (p  -  j)  -f-,  equals  p  and  that 
i'  -  j  equals  (»'  -  j  -  1)  -  1.  equation  A-7  can  be  rewritten  as: 


=  ^  \  |  pB  i  /  -  j  - 1.  p  -  j) 

/  t  \  r(if  -  j  -  i)f(p  -  j) 

\  j  I  pr  (i*  -  j  -  i)  -  (p  ~  j) 

(<')'■  (t'-jVT(p-j) 

V'-jYJ 'P  -p)  ■ 


B{  )  and  T(  )  are  the  standard  Beta  and  Gamma  functions 

Substituting  the  RHS  of  equation  A-8  for  p ,,  in  equation  in  of  the  main  text  yield 
the  following  result  for  the  exponential  case: 


Q  -  En=.  ,  Cn  ~ 


(  m  -  n  ]'  1  T  ■ ,  -  m-  k  ] 
im~k)'  -pr 

,  -  1  m  -  n; '  *  1 


( m  -  k  | '  P  f ;  m  -  k~  1 ) 


.  k  =  0.1.  ..m  -  1 


where  (m  -  n)’  =  min{m  -  n  —  l.m}. 

(.'sing  the  fact  F(p  ~  rn  -  k  —  I)  =  (p~  m  -  k)F(p  -  m  -  k)  and  cancelling  out  -~-k 
and  r  '  ■  equation  A-9  may  be  rewritten  as: 


(A- 10 


^  _  \^k  m  -  n  ■ '  '  F  I,.  -  m  -  k  i 

*  «— 1 =  n  <  rr,  -  k  P  f  ,  - 


:  k  -  0.1. 
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For  k  =  0.  equation  A- 10  yields: 


C  -  C.  %pr~,  ,~rr‘— 

m  '  f , ,  -  rr  *  1 


p  —  rn  -  0 


rn  -  0 


Combining  equations  A- 1C'  and  A- 11  fields 


c 

=  1. 

(A  -12 

m  -  0 

Cl 

=  — — — C.  . 

P 

Clc-l 

_  * 

(  -  rr<  -  k 


k  —  1.2 . m  -  1 

The  summation  of  expressions  involving  F  p  -  (m  —  n)1  —  1  is  analytically  difficult 
therefore,  the  following  approach  is  taken: 


m  -  (A*  -  1)  '  =  (m  -  k)  -  1  '  =  (m  -  k  -  1  )(m  ~  At)! 


Therefore 


(m  A-)’  = 


m  -  [k  -  1  l  ' 
m  -  k  -  1 


(A  13 


similarly 


Tp  -  m  -  (1-  -  1)  =  T  ( p  -  m  -  A* ;  -  1 

=  ( p  ~  m  -  A- )  f  ( p  ~  n>  -  A- ) 


Therefore: 


Y(p  -  m  -  A-)  = 


T  p  ~  m  (A-  1 ) 

p  —  ni  -  k 


(A  11 


l  sing  equations  A-13  and  A- 1 4 .  the  summation  component  of  equation  A-13  can  be 


rew  rit  ten  as 


Scrutinizing  the  RHS  of  equation  A-lo  with  the  help  of  equation  A-*  Yields: 


('  (m  ~  rt)>  '  f  {p  -  m  -  k) 

^  (”>  -  k)'  (m  -  n)'  -  1  -  p 


m  -  k  ~  1 
p  —  m  -  k 


k 


m  - '  k  -  1  i 


(A- 16) 


In  order  to  transform  the  summation  component  in  the  RHS  of  equation  A-16  into 
a  more  useful  form,  let  us  rewrite  equation  ID  of  the  main  text  as: 


C 

Cr 


c 


k  -  1 


_  V*-  1  (' 

—  n~  .  f. 


P'n  -  r  r»; 


Pm  -  k  m  -  *  -  1 


./•  =  1.2. 


(A- 17) 


Multiplying  both  sides  by  adding  n  -  k  to  the  summation  range,  and 

moving  the  summation  to  the  other  side,  equation  A-17  yields: 
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l  tilizing  equations  A-16  and  A- 1 1? .  equation  A-I2  can  be  rewritten  as 


a 


a 
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fi 
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If  m  =  0.  the  system  is  of  no  interest  at  all.  If  m  =  1.  the  system  is  completely 
described  by  C  and  C] .  If  m  >  2.  equation  A-19  with  k  =  1  reduces  to: 


m  -  0 


C, 
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Again,  if  m  =  2.  the  system  is  completely  described  by  CV  .  Cj.  and  C;.  If  m  >  3 
equation  A-19  with  fc  =  2  reduces  to: 


m  0 

- ( 

P 

m  -  1 

- < 
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/->  m  -  2  —  1  /'•’ 
^  2  “  - 1 


r  -  m  -  » 
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v-m-  2 


p  —  m  -  2  „  ni  —  1 
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p  p 


ic,  -  C,  =  —  C; 


By  now.  the  following  recursive  relationship  is  suspected: 


C-,  =  1 


(A-20; 
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Assume  that  equation  A-20  holds  for  k  =  0.1 . A-  -  1  for  an  arbitrary  A-'  To 

complete  the  induction  process,  it  must  be  shown  that  equation  A-20  holds  for  A*  —  A-' 
If  A'  =  m.  the  process  is  complete.  If  A-’  £  m  -  1.  equation  A-19  with  A-  -  k'  yields: 


CV_, 


cv 


m  -  k"  -  1 
r-m-  k‘  -  l 


i—m-k' 


rn 


l 


(A-21 ) 


Because  it  is  assumed  that  equation  A-20  holds  for  k  -  k'  -  1.  that  is.  Ck-  = 
— -Cf-i-  equation  A-21  can  be  rewritten  as; 


p  —  m  -  A-'  „  m-A*' 

C*.-1  =  - cv  -  c*.  =  - Ck- 


which  completes  the  induction  process. 


Having  established  the  validity  of  the  relationship  given  by  equation  A-20  for  all  As 
in  the  appropriate  range,  it  should  be  noted  that  an  equivalent  form  is: 


or: 


C.  -  I 


( A — 2  2 ) 


cv,  =- 


m  -  1 
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C\.  =  1 


Ck  = 


(m  -  A)'/ 


:  A*  —  1.2. 


1 
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Substituting  the  RHS  of  equation  A-23  for  the  Cs  in  equation  9  of  the  main  text 
yields  the  following  steady-state  probabilities  for  the  exponential  case: 


m  —  ■  m-£)  ' 
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—  *  :  m  —  n ;  Vr  — ' 


A-  =  0.1. 


y  m  j 
— *  n  =  • ‘  it 


(A -24) 


Cancelling  out  and  ^ .  looking  at  (m  -  n)  as  one  index,  and  reversing  its  order 
of  summation  yields: 


pk  k'. 

r*  =  m 


which  are  the  steady-state  probabilities  for  an  A/  A/  rn  m  o c  queueing  system. 
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